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SECTION  I 
PRELUDE 


INTRODUCTION 

This  report  documents  a  study  concerned  with  the  use  of  a  centrifuge  as 
an  experimental  device  on  which  free-field  blast  parameters  could  be  evalu¬ 
ated.  The  centrifugal -simulation  of  high-explosive  (HE)  charges  was  investi¬ 
gated  as  a  means  of  improving  the  present  techniques  used  to  predict  free- 
field  blast  parameters.  The  blast  parameters  of  interest  included  normal 
stress,  particle  velocity,  and  acceleration.*  However,  because  of  a  lack  of 
suitable  accelerometers  and  velocity  gages,  this  study  is  concerned  only  with 
normal  stress. 

The  gravity  field  created  by  the  rotating  centrifuge  increase  the  effec¬ 
tive  body  force  of  gram-weight  charges,  thus  simulating  a  larger  explosive. 

The  prototype  explosive  is  a  cylindrical-shaped  ,  cased  explosive  embedded  in 
a  concrete  burster  slab.  The  prototype  burster  slab  is  supported  by  a  deep 
stratum  of  sand.  A  definition  of  the  variation  of  normal  stress  in  respect  to 
distance  in  the  region  near  the  explosive  was  a  primary  goal  of  this  study. 

BACKGROUND 

Recent  developments  in  precision-guided,  conventional  weapons  have 
greatly  increased  the  vulnerability  of  military  facilities.  It  is  now  pos¬ 
sible  for  an  attacker  to  precisely  deliver  relativly  large  conventional  weap¬ 
ons  on  selected  targets.  One  means  of  protecting  strategic  facilities  from 
such  weapons  is  a  concrete  burster  slab  and  intervening  sand  stratum.  Not 
only  can  weapons  be  accurately  directed  at  these  sites,  but  contemporary 
weapons  can  penetrate  part  of  the  burster  slab  before  detonating.  Economic 


*In  this  report  the  term  pressure  refers  to  the  detonatioh  pressure  of  an 
explosive.  The  pressure  in  the  free  field  as  a  result  of  detonating  the 
weapon  is  referred  to  as  stress.  No  distinction  is  made  between  hydrodynamic 
or  soil  stress. 


design  and  survival  of  buried. structures  requires  a  thorough  knowledge  of  the 
stress  loadings  created  by  the  weapon  when  it  detonates.  The  assumption  is 
made  that  the  burster  slab  defeats  the  weapon;  that  is,  the  weapon  lodges  in 
the  burster. slab  and  then  detonates. 

Extensive  theoretical  and  experimental  research  was  performed  after  World 
War  II  (Reference  1)  to  describe  HE-generated  blast  waves  and  their  propaga¬ 
tion  and  attenuation  through  various  media.  However,  much  of  that  work 
involved  the  use  of  spherical,  uncased,  surface-detonated  charg  Moreover, 
most  of  the  data  associated  with  those  early  tests  were  collect  at  rela¬ 
tively  large-scaled  distances  (X)  because  of  the  instrunentatio  difficulties 
associated  with  the  near  field,  X  <  2.  Most  of  the  work  done  d  the  1950s 
and  1960s  was  in  the  field  of  nuclear  (or  HE  simulation  of  nuc'i  j  olast 
environments.  Virtually  none  of  these  programs  included  a  burster  slab,  and 
only  a  few  of  the  tests  included  cased  explosives.  Model  tests  (circa  1980) 
have  indicated  that  cased  explosives  which  were  detonated  within  a  burster 
slab  did  not  yield  normal  stress  attenuation  rates  consistent  with  those 
extrapolated  from  the ^vffies^maTticined  afiove.  ^7 

PURPOSE 

The  purpose  of  the  effort  documented  in  this  report  was  to  develop  an 
experimental  technique  to  define  the  free-field  blast  parameters  in  the  region 
near  cylindrical-shaped,  cased,  HE  charges,  in  soil  and  concrete-soil  media. 

SCOPE 

To  conduct  this  research  it  was  necessary  to  (1)  review  the  current  state 
of  the  art  regarding  the  definition  of  blast  wave  generation  and  propagation 
through  soil  and  concrete-soil  media,  (2)  define  the  parameters  necessary  to 
describe  such  phenomena,  (3)  review  the  experimental  techniques  currently 
employed  for  such  studies,  (4)  examine  the  scaling  laws  related  to  blast  wave 
propagation,  and  (5)  critically  evaluate  the  instrumentation  and  data  acquisi¬ 
tion  systems  used  in  these  types  of  studies.  Subsequent  to  these  reviews,  the 
experimental  technique  recommended  to  meet  the  purpose  of  this  study  was  to  be 
evaluated. 


CRITERIA 


The  experimental  technique  recommended  to  satisfy  the  objectives  of  this 
study  was  required  to  have  the  following  capabilities: 

1.  Must  allow  for  variation  of  charge  weight.  The  maximum  scaled  charge 
could  not  exceed  15  pounds  of  trinitrotoluene  (TNT)  to  simulate  a  full-scale 
TNT  charge  weight  of  1000  pounds. 

2.  Must  allow  variation  of  the  depth  of  burial  ud  to  a  scaled  distance 
of  3  feet. 

3.  Must  allow  variation  of  the  charge  orientation  from  the  normal  to  an 
angle  of  60  degrees. 

4.  Must  allow  variation  of  the  scaled  thickness  of  the  prototype  con¬ 
crete  burster  slab  ranging  from  3  to  6  feet. 

5.  Must  allow  representation  of  the  concrete-soi 1  media  as  a  half 
space. 

6.  Must  provide  quantifiable  data  on  the  blast  parameters  in  the  scaled 
distance  range  of  0.5  to  2  in  the  soil  and  concrete  soil  media. 

7.  Must  be  repeatable  or  reusable  with  minor  refurbishing. 

APPROACH 

The  initial  effort  in  this  research  was  principally  concerned  with  a 
critical  review  of  blast  parameter  evaluation  techniques  with  respect  to  their 
potential  to  meet  the  full  range  of  capabilities  outlined  above.  During  a 
briefing  held  at  the  Air  Force  Engineering  and  Services  Center  on  1  April 
1981,  the  author  outlined  several  potential  experimental  techniques  that  could 
be  used  as  simulators.  These  included  a  modified  Stanford  Research  Institute 
(SRI)  single-source  line  loader,  small-charge  field  tests,  very  small-charge 
model  tests,  and  the  use  of  a  centrifuge.  Of  these  techniques,  only  the  cen¬ 
trifugal  concept  had  the  potential  to  meet  all  of  the  preceding  criteria. 

Following  that  briefing,  the  scope  of  this  project  was  extended  to 
include  an  experimental  evaluation  of  the  centrifugal  concept.  This  report 
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documents  that  study.  Section  II  of  this  report  is  a  discussion  of  the  appro 
priate  scaling  laws  derived  for  gravity  scaling.  All  aspects  of  the  experi¬ 
mental  phase--pressure  gage  development,  data  acquisition,  model  explosives, 
system  modeling,  and  an  outline  of  the  experimental  program— are  discussed  in 
Section  III.  Experimental  results  are  presented  in  Section  IV.  Section  V  is 
a  discussion  of  those  results.  An  extensive  bibliography  is  included  at  the 
end  of  this  report. 


SECTION  II 

DIMENSIONAL  ANALYSIS  AND  SIMILITUDE 


DIMENSIONAL  ANALYSIS 

The  techniques  of  dimensional  analyses  are  used  in  experimental  programs 
in  which  the  difficulty  encountered  in  developing  a  closed-form  solution  (dif¬ 
ferential  equation)  which  describes  the  phemonena  is  so  complex  as  not  to  be 
readily  discernible.  At  best,  such  equations  (when  they  can  be  developed)  are 
based  upon  first-order  concepts.  Thus  parametric  studies  based  upon  differen¬ 
tial  equations  suggest  significantly  different  system  behavior  than  that 
observed  experimentally.  This  departure  is  due  to  the  imperfection  of  human 
observations  (objective  uncertainty),  the  imperfection  of  intellectual  con¬ 
cepts  devised  to  reproduce  physical  phenomena  (subjective  ignorance),  unreal¬ 
istic  mathematical  assumptions  (used  to  ease  the  computational  effort),  and 
experimental  error  in  the  test  program.  Error  in  the  test  program  exists 
regardless  of  the  analyses  employed.  However,  it  is  more  difficult  to 
identify  error  when  the  data  are  interpreted  using  first-order  concepts. 

Subjective  ignorance  disappears  in  the  formulation  of  a  problem  based 
upon  the  concepts  of  dimensional  analysis.  This  is  because  the  functional 
form  of  the  relationship  between  the  dependent  and  independent  variables  is 
determined  experimentally.  Thus  the  weakness  of  dimensional  analysis  is  imme¬ 
diately  obvious;  a  class  or  family  of  solutions  (parametric  studies)  cannot  be 
generated  as  in  the  case  of  closed-form  solutions.  Rather,  each  functional 
relationship  must  be  determined  through  experimentation. 

Objective  uncertainty  also  exists  in  dimensional  analysis,  but  the  incom¬ 
pleteness  that  such  uncertainty  suggests  is  immediately  obvious  in  an  examina¬ 
tion  of  the  experimental  data.  Thus,  once  the  possibility  of  significant 
constant  or  systematically  changed  errors  is  identified  and  removed  from  the 
experiment  or  minimized,  and  the  remaining  errors  are  primarily  random,  the 
deficiencies  in  the  model  due  to  objective  uncertainty  can  be  identified. 

By  the  simple  process  of  considering  additional  independent  variables  the 
model  can  be  revised,  as  often  as  necessary,  until  a  satisfactory  correlation 
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is  achieved  between  the  dependent  variable  and  the  independent  variables.  .The 
process  is  actually  one  of  eliminating  objective  uncertainty  from  the  analy¬ 
sis.  In  a  sense,  dimensional  analysis  contains  a  feedback  mechanism  which 
permits  the  analyst  to  eventually  identify  the  significant  variables  in  the 
problem  under  consideration.  In  short,  the  technique  greatly  minimizes  objec¬ 
tive  uncertainty.  The  converse  is  the  situation  in  closed-form  solutions.  In 
these  cases,  experimental  data  always  depart  from  the  mathematical  model 
because  of  random  errors  in  the  test  procedures  and,  more  importantly,  because 
the  functional  relationship  was  established  a  priori .  In  essence,  the  researcher 
uses,  an  Incomplete  theory  to  describe  a  physical  event.  The  data  from 
the  physical  event  are  manipulated  to  conform  to  the  incomplete  concept  (objec¬ 
tive  uncertainty),  and,  in  so  doing,  the  behavior  of  the  event  is  modified  by 
subjective  ignorance.  The  physical  system  is  forced  to  respond  in  a  manner 
that  is  not  compatible  with  its  essential  boundary  conditions  or  intensive 
material  properties.  The  great  tragedy  is  that  the  results  of  these  force 
fits  often  become  engineering  standards.  Hopefully,  none  of  the  existing 
strategic  facilities  designed  on  the  basis  of  these  standards  will  ever  be  put 
to  the  test. 

Figure  1  is  an  orthogonal  view  of  the  prototype  burster-slab  sand  system. 
Also  shown  is  the  centrifugal  model  of  the  system.  An  immediately  obvious 
problem  with  the  model  is  that  its  boundaries  provide  a  surface  for  reflection 
of  the  explosive-induced  pressure  wave.  The  size  of  the  model  is  limited  to 
space  provided  on  the  centrifuge  test  platform.  Because  the  boundary  is  nec¬ 
essary,  data  traces  will  need  to  be  examined  for  boundary  effects. 

The  primary  dependent  parameter  of  interest  in  this  study  is  normal 
stress  o,  (ML"1!*2).  Other  possible  dependent  parameters  are  acceleration, 
particle  velocity,  displacement,  or  characteristic  time.  A  list  of  the 
assumed  (objective  uncertainty)  independent  parameters  is  presented  in 
Table  1. 
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Sand 

Pressure  sensors 
Centrifuge  base 


Burster  slab 


TABLE  1.  INDEPENDENT  PARAMETERS  FOR  THE 

CENTRIFUGAL  MODELING  OF  BLAST  PRESSURE. 


Parameter 

Description 

Dimensions 

Explosives 

Energy,  EQ 

ml2t-2 

Pressure,  PQ 

ML-1T-2 

Depth  of  Burial,  D 

L 

Burster-Si  ab 

Thickness,  H 

L 

Mass  Density,  px 

ML"3 

Dilatational  Wave  Speed,  Cx 

LT-1 

Poisson's  Ratio, 

— 

Strength, 

ML-1T~2 

Soil  (Sand) 

Mass  Density,  p2 

ML-3 

Dilatational  Wave  Speed,  C2 

LT-1 

Poisson's  Ratio,  v2 

— 

Strength  Parameter  a2 

ML-1T"2 

Other 

Gravity,  G 

LT-2 

Parameters 

Bucket  Dimensions,  r,d 

L 

Range 

L 

The  Buckingham  ^-theorem  (Reference  2)  is  fundamental  to  the  concept  of 
dimensional  analysis.  The  "-theorem  is  adequately  presented  in  countless 
books  concerned  with  hydraulics  and  fluid  mechanics;  its  implementation  is  not 
presented  here,  only  the  results.  An  excellent  treatment  of  the  "-theorem 
applied  to  ground  motion  problems  is  included  in  Reference  3.  It  is  suffi¬ 
cient  to  state  that  by  means  of  the  w -theorem  the  independent  parameters  are 
formed  into  ratios,  called  "-terms.  A  significant  aspect  of  each  "-term  is 
that  it  is  a  nondimensional  term.  Each  term  is  an  independent  variable  ratio 
which,  when  combined  with  all  other  "-terms,  defines  the  functional  relation¬ 
ship  between  the  dependent  and  independent  variables.  Thus,  nondimensional 
stress  is  expressed  as 


where  ?1(  or  c2»  is  any  characteristic  length,  such  as  radial  position, 
burster-slab  thickness,  any  dimension  of  the  explosive,  or  any  other  charac¬ 
teristic  length.  The  dots  indicate  that  Equation  (1)  is  not  complete.  Com¬ 
pleteness  might  be  achieved  by  the  addition  of  * -terms  made  up  of  such  vari¬ 
ables  as  sand  moisture,  viscosity,  or  strain-rate  effects.  The  concept  of 
objective  uncertainty  is  readily  admitted  to  this  point.  Only  an  experimental 
exercise  will  produce  the  data  to  define  the  functional  form  of  Equation  (1); 
and  an  analysis  of  variance  will  suggest  the  need  for  additional  terms  in  the 
equation. 

SIMILITUDE 

For  data  from  model  tests  to  be  useful  in  predicting  prototype  behavior, 
'invariance  must  exist  between  the  model  and  the  prototype.  Because  a  model  is 
used,  it  becomes  obvious  that  its  dimensions  must  be  proportional  (similar)  to 
those  of  the  prototype.  The  scaling  relationships  used  to  provide  similarity 
between  the  model  and  prototype  are  derived  from  dimensional  analysis.  Simil¬ 
itude  involves  scaling  not  only  characteristic  lengths  (dimensions)  but  also 
time,  material  properties,  and,  in  the  case  of  this  study,  gravity.  Thus  any 
characteristic  length  related  to  the  prototype  is  divided  by  n  to  yield  the 
corresponding  dimension  for  the  model.*  The  scaling  relationships  derived  for 
this  study  are  presented  in  Table  2. 

These  relationships  are  derived  from  the  * -terms  developed  through  appli¬ 
cation  of  the  concepts  of  dimensional  analysis.  Some  of  these  u-terms  are 
included  here  to  illustrate  the  concept. 

q 

*i  =  p 2^2  Equivalency  of  stress 


"2  -M& 


IT,  * 


1/3 


Geometric  similarity  of  charateristic  length 


Similarity  of  medium  properties 


*In  this  report  n  is  the  ratio  of  the  centrifugal  acceleration  to  the  gravita¬ 
tional  field  at  the  earth's  surface.  For  simplicity  this  ratio  is  called  the 
G-level,  implying  n  =  10  if  the  centrifugal  acceleration  is  322  ft/s2. 


TABLE  2.  SCALING  RELATIONSHIPS  FOR  GRAVITY  SCALING. 


Quantity 

Full  Scale 

Model  Scale 

Linear  Dimension 

1 

■SM 

Area 

1 

WmSmm 

Volume 

1 

Time 

1 

l/n 

Velocity 

.  1 

1 

Acceleration 

1 

n 

Mass 

1 

l/n3 

Force 

1 

l/n2 

Energy 

1 

l/n3 

Stress 

1 

1 

Strain 

1 

1 

Density 

1 

1 

Gravity  scaling  of  blast  energy 


Gravity  scaling  of  time 


Eighteen  ir-terms  are  associated  with  the  independent  parameters  listed  in 
Table  1  and  the  five  response  (dependent)  parameters.  The  scaling  factors 
presented  in  Table  2  are  developed  from  these  "-terms. 


Considering  gravity  scaling  of  energy  for  example: 


Therefore,  energy  scales  as  1/n3 . 

The  "-term  of  importance  to  this  project  is  "4.  This  term  can  be  written 
in  terms  of  mass  as 


Q  *  heat  of  detonation  per  unit  mass  of  explosive 
<S  =  initial  density  of  the  explosive 
W  *  mass  of  the  explosive 
G  *  gravity 

By  equating  the  *4-term  of  the  model  to  that  of  the  prototype,  the 
similitude  requirement  for  explosives  can  be  developed.  This  development  is 
discussed  in  Section  III,  Explosives.  Further  details  on  modeling  are  also 
presented  in  Section  III,  Experimental  Program. 
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SECTION  III 

EXPERIMENTAL  PROCEDURES 


INTRODUCTION 

The  experimental  phase  of  this  study  consisted  of  detonating  a  small 
explosive  charge  embedded  in  a  model  burster  slab  and  recording  the  free-field 
normal  stress  at  selected  locations  in  a  contiguous  sand  medium.  The  burster 
slab  and  sand  continuum  were  contained  in  a  22-inch  diameter  by  a  12-inch  high 
aluminum  test  fixture  which  was  mounted  on  the  test  platform  on  a  30,000  G-lb 
centrifuge.  The  explosive  charge  was  detonated  with  the  centrifuge  rotating 
at  a  predetermined  constant  angular  velocity.  The  gravity  field  which  devel¬ 
ops  because  of  the  rotational  effects  of  the  centrifuge  increases  the  effec¬ 
tive  body  force  of  gram-weight  charges,  thus  simulating  a  larger  explosive. 
This  technique  permits  the  yield  of  a  small  explosive  to  simulate  that  of  a 
large  explosive.  The  benefits  of  such  an  approach  are  obvious.  A  series  of 
22  of  these  experiments  was  performed  as  described  later  in  this  section. 

To  perform  these  experiments  it  was  necessary  to  (1)  become  familiar  with 
the  centrifuge  and  its  operation,  (2)  develop  safety  procedures  relative  to 
its  operation  and  the  detonation  of  the  explosives,  (3)  establish  the  geometry 
and  weight  of  each  explosive,  (4)  model  the  burster  slab  and  sand  media,  and 
(5)  develop  a  suitable  stress  gage  and  associated  data  acquisition  system. 

CENTRIFUGE 

The  centrifuge  used  to  conduct  the  experimental  phase  of  this  project  is 
a  Genisco  model  E185  (Serial  Number  11).  A  photograph  of  this  centrifuge  is 
presented  in  Figure  2.  The  boom  is  nominally  6  feet  long.  The  centrifuge  is 
situated  below  floor  level  in  a  concrete  pit  located  in  Building  1001  on  Kirt- 
land  Air  Force  Base.  The  centrifuge  is  covered  with  a  reinforced  concrete 
roof  with  access  to  the  pit  through  steelplate  doors.  The  centrifuge  is  con¬ 
trolled  by  suitable  electrical  and  hydraulic  controls  located  aboveground  in  a 
console  remote  to  the  centrifuge.  The  centrifuge  was  originally  built  to  test 
aviation  electronic  and  mechanical  systems  under  G-loadings  similar  to  those 
experienced  during  flight.  A  variable-speed  hydraulic  drive  system  is  used  to 
rotate  the  centrifuge.  The  hydraulic  power  unit  is  enclosed  in  the  console; 


the  fluid  motor  is  installed  at  the  bottom  of  the  rotor  assembly.  The  boom 
consists  of  two  symmetrical  cantilever  arms,  each  of  which  supports  a  test 
platform.  The  platforms  are  30  inches  square  and  were  designed  to  operate  in 
one  of  five  fixed  positions:  horizontal,  45,  90,  135,  or  180  degrees 
(inverted).  Electrical  interlock  switches  prevent  rotation  of  the  centrifuge 
unless  the  platforms  are  securely  locked  into  one  of  the  specified  positions. 
However,  for  this  research  the  platforms  were  modified  to  permit  the  centri¬ 
fuge  to  operate  in  the  swing-bucket  mode  rather  than  fixed-bucket  as 
described.  These  modifications  were  made  in  accordance  with  the  instructions 
contained  in  a  report  prepared  by  the  manufacturer  for  the  New  Mexico  Engi¬ 
neering  Research  Institute  (NMERI)  (Reference  4).  The  swing-bucket  modifica¬ 
tion  was  necessary  to  facilitate  handling  the  test  fixture  which  weighed,  on 
the  average,  365  pounds. 

In  the  swing-bucket  mode,  the  test  fixture  could  be  easily  lowered  onto 
the  horizontal  test  platform  and  then  permitted  to  swing  into  the  90-degree 
position  as  the  centrifuge  was  subsequently  accelerated  to  a  prescribed  con¬ 
stant  angular  velocity.  Mechanical  scratch  gages,  mounted  on  the  boom  and 
installed  to  make  contact  with  one  of  the  arms  of  the  swinging  platforms, 
indicated  that  the  90-degree  position  was  consistently  achieved. 

The  centrifuge  is  equipped  with  28  slip  rings  which  provide  for  elec¬ 
trical  connections  between  the  rotating  test  object  and  the  remote  control 
console.  The  boom  is  free  to  pivot  approximately  1  degree  about  a  horizontal 
transverse  axis.  Static  balancing  is  achieved  if  identical  test  objects  are 
placed  on  each  of  the  opposite  test  platforms.  If  a  single  object  is  to  be 
tested  (as  in  this  research)  counterweights  are  mounted  on  the  opposite  test 
platform.  An  automatic  dynamic  balancing  feature  of  the  machine  eliminates 
the  necessity  for  static  balancing  to  less  than  10  pounds.  Wow  and  drift  of 
the  boom  assembly  during  a  1-minute  period  is  less  than  0.5  percent  of  the 
operating  rate  at  any  speed.  Speed  is  indicated  on  a  direct  reading  magnetic- 
impulse-type  tachometer  indicator,  accurate  to  within  1  1/4  percent  of  full 
scale.  Precision  speed  indicators  are  provided  by  a  pulse  generator  unit  and 
an  electronic  counter  on  a  four-decade,  digital  display.  Accuracy  of  the  unit 
is  within  0.1  percent  at  angular  velocities  above  12  rpn.  The  average  angu¬ 
lar  speed  over  a  1-  or  10-second  interval  is  indicated  directly  in  revolutions 
per  minute. 
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Two  test  objects„each  weighing  up  to  500  pounds,can  be  tested  simultan¬ 
eously,  one  on  each  of  the  two  test  platforms.  The  maximum  range  of  the  cen¬ 
trifuge  is  100  G  at  a  72-inch  radius,  22  to  220  rpm.  infinitely  variable. 

The  centrifuge  is  rated  at  30,000  G-lb.  Thus,  at  its  maximum  centrifugal 
acceleration  (100  G)  a  load  of  300  pounds  can  be  applied  to  each  test  plat¬ 
form.  An  acceleration  of  60  G  would  be  the  maximum  allowed  if  500  pounds  were 
applied  to  each  test  platform. 

SAFETY  PROCEDURES 

The  centrifuge  is  placed  in  a  circular  reinforced  concrete  pit  approxi¬ 
mately  7  feet  deep.  The  roof  of  the  pit  is  also  made  of  reinforced  concrete. 
Several  steel  pi  ate  doors  mounted  in  the  roof  provide  access  to  the  floor  of 
the  pit.  Thus,  personnel  are  adequately  protected  from  the  inadvertent  emis¬ 
sion  of  missiles  from  the  rotating  centrifuge.  Interlock  switches  on  the 
access  doors  serve  as  safety  devices.  These  prevent  the  rotor  assembly  from 
operating  when  any  door  is  not  properly  seated  against  its  jam.  Conversely, 
these  switches  initiate  an  immediate  shutdown  of  the  rotor  if  an  attempt  is 
made  to  open  a  door  while  the  centrifuge  is  rotating.  Other  similar  devices 
prevent  the  rotor  from  exceeding  a  minimum  angular  velocity  if  a  dynamic 
unbalance  exists,  or  if  the  weight  of  the  test  object  exceeds  500  pounds. 
Further  details  concerning  the  Genisco  centrifuge  and  its  operation  or  main¬ 
tenance  are  contained  in  Reference  5. 

The  placement  of  the  test  fixture  containing  the  burster  slab,  with  the 
stress  gages  embedded  in  the  sand,  was  under  the  direction  of  the  project 
engineer.  However,  once  the  fixture  was  secured  to  the  platform,  an  ordnance 
technician  served  as  the  test  conductor.  This  person  took  charge  of  the  oper¬ 
ation  and  served  as  the  responsible  safety  officer  during  the  arming  and  deto¬ 
nation  of  the  explosive.  This  also  included  a  postshot  inspection  of  the 
centrifuge  pit  to  ensure  that  the  area  was  free  of  all  explosive  hazards.  An 
Explosives  Facility  License,  Reference  6,  was  granted  to  the  Civil  Engineering 
Research  Facility  to  detonated  gram-weight  explosives  in  the  centrifuge  facil¬ 
ity.  The  approved  Centrifuge  Countdown  Procedure  is  presented  as  Appendix  A. 


EXPLOSIVES 


In  simulation  testing,  an  attempt  is  made  to  maintain  similitude  between 
the  full-scale  object  (prototype)  and  a  scaled  model  of  the  object.  The  con¬ 
cept  evaluated  in  this  project  is  ba.-ed  on  the  consideration  that  a  small 
explosive  charge  (model)  subjected  to  an  acceleration  field  can,  upon  detona¬ 
tion,  simulate  the  effect  of  a  larger  explosive  (prototype).  The  centrifuge 
is  used  to  create  the  acceleration  field.  In  regard  to  the  prototype  explo¬ 
sive,  simulation  requires  that  the  characteristic  dimensions  and  yield  of  the 
model  be  determined  according  to  the  appropriate  scaling  relationships  (Sec¬ 
tion  II).  These  relationships  require  that  the  characteristic  lengths  of  the 
model  scale  as  the  reciprocal  of  the  G-level;  and  that  the  energy  (yield) 
scales  as  the  reciprocal  of  the  G-level  cubed,  with  respect  to  those  values 
for  the  prototype.  Thus,  if  a  particular  weapon  (prototype)  is  to  be  simula¬ 
ted,  the  application  of  these  scaling  relationships  yields  the  necessary 
dimensions  and  explosive  weight  for  the  model.  It  is  obvious  that  an  infinite 
number  of  models  exist  for  a  particular  weapon  if  the  G-level  is  infinitely 
variable  as  with  a  centrifuge.  This  is  advantageous  because  it  permits  the 
consideration  of  a  range  of  models  in  simulation  experiments. 

This  flexibility  is  also  significant  in  that  it  improves  the  outlook  for 
full  experimental  similitude  when  other  secondary  test  parameters  are  consid¬ 
ered.  For  example,  the  scaling  relationships  may  suggest  that  the  size  of 
other  features  included  in,  or  necessary  to,  the  experiment  be  so  small  as  not 
to  be  commercially  available,  difficult  to  manufacture  or  handle,  or  simply 
not  available.  Thus,  to  perform  the  intended  experiment,  the  G-level  is 
established  on  the  basis  of  available,  achievable  model  dimensions,  or  those 
of  the  secondary  features  of  the  prototype  system.  The  advantage  of  an  infi¬ 
nitely  variable  G-field  (centrifuge)  is  that  it  provides  a  means  of  approach¬ 
ing  the  ideal  condition  of  constructing  an  experiment  that  is  in  compliance 
with  all  conditions  required  for  similitude.  However,  this  condition  is 
rarely  met,  not  only  because  of  physical  limitations,  but  because  of  the 
analyticity  from  which  the  scaling  laws  are  established.  By  their  very  nature 
these  laws  are  subjective  in  that  they  are  derived  through  an  interrogative 
process  (dimensional  analysis)  that  is  highly  subjective  and  admittedly  incom¬ 
plete.  The  process  (dimensional  analysis)  is  used  when  science  admits  to  its 
inability  to  develop  an  exact  mathematical  statement  for  the  problem.  These 


limitations  notwithstanding,  it  is  not  suggested  that  the  use  of  dimensional 
analysis  be  abandoned.  To  the  contrary,  a  correctly  formulated  problem  pro¬ 
vides  insight  concerning  the  parametric  relationships  which  exist  between  the 
dependent  and  perceived  independent  variables.  These  relationships  serve  to 
improve  the  understanding  of  the  problem  and  to  suggest  experimental  and  ana¬ 
lytical  tactics  which  rarely  are  achieved  through  physical  reasoning  alone. 
Perfect  correlations  cannot,  however,  be  expected  in  experimental  observations 
when  the  results  are  formulated  in  terms  of  dimensional  parameters,  because 
the  theory  is  always  incomplete,  because  full  simulation  cannot  be  achieved, 
and  because  of  experimental  errors. 

The  limitations  and  considerations  discussed  above  are  fully  demonstrated 
in  the  details  leading  to  the  selection  of  model  explosives  for  use  in  this 
project.  While  a  particular  weapon  (prototype)  was  not  identified  with  this 
project,  the  weapon  considered  was  identified  as  a  general-purpose  bomb  (GPB) 
in  the  total  weight  range  of  500  to  2000  pounds.  On  the  surface,  simulation 
appears  easy;  application  of  the  appropriate  scaling  laws  to  the  explosive 
weight  and  dimensions  of  these  bombs  yields  the  corresponding  values  for  the 
model  explosive,  the  G-level  known  or  selected  a  priori .  For  example,  a  2000- 
pound  GPB  containing  1117  pounds  of  TNT  could  be  simulated  by  107,  13,  4,  or 
1.7  grams  of  <^lotriniet^yl^etrinitramine  (RDX)  when  accelerated  to  20,  40, 

60,  or  80  G  respectively.  However,  the  actual  weight  of  TNT  in  any  GPB  is 
variable,  depending  upon  the  manufacturing  processes.  These  processes  affect 
not  only  the  properties  of  the  TNT,  but  also  influence  the  shape  of  the  explo¬ 
sive  contained  within  the  bomb  case.  Likewise,  the  dimensions  and  shape  of 
the  case  varies,  as  well  as  its  intensive  properties. 

Complete  simulation  requires  a  review  of  all  of  the  above  factors.  Sta¬ 
tistical  uncertainties  are  obviously  present  in  respect  to  both  the  intensive 
and  extensive  properties  of  the  prototype  to  be  simulated.  True  dimensional 
and  energy  scaling  would  lead  to  model  charges  that  would  be  one  of  a  kind  and 
expensive  to  fabricate.  Not  only  would  the  cost  of  the  project  be  signifi¬ 
cantly  affected,  but  it  is  doubtful  that  such  detail,  though  desirable,  would 
enhance  the  real  purpose  of  this  program,  namely,  to  demonstrate  the  feasibi¬ 
lity  of  the  centrifuge  as  a  simulator  to  define  blast  parameters  in  the  region 
near  a  high-explosive  charge  placed  in  a  concrete-soil  medium. 


An  alternate  approach  was  taken,  though  it  was  still  guided  by  the  appro¬ 
priate  scaling  laws.  This  consisted  of  dealing  with  what  could  be  considered 
the  average  properties  of  general-purpose  bombs.  The  first  consideration  was 
that  of  determining  the  weight  of  a  model  RDX  charge,  expressed  as  a  function 
of  the  G-level,  required  to  simulate  a  particular  general-purpose  bomb.  Equa¬ 
tion  (2)  was  used  for  this  purpose.  The  results  of  these  computations  are 
shown  in  Table  3. 


TABLE  3.  THEORETICAL  MODEL  EXPLOSIVE  SIMULATION  WEIGHTS. 


Bomb 

Designation, 

pound 

60  g 

80  g 

100  g 

Weight  of  RDX  in  grams 

250 

12.35 

1.54 

1 

m 

B9 

■ 

25.57 

3.20 

EH 

EH 

53.43 

6.68 

1.98 

m 

1 

106.95 

13.37 

3.98 

1.67 

HI 

These  results,  based  upon  similitude  concepts,  indicate  that  small 
charges,  when  subjected  to  the  indicated  acceleration,  can  simulate  the  yield 
of  large  weapons.  Because  of  the  uncertainties  related  to  centrifuge  explo¬ 
sive  testing,  the  NMERI  Safety  Officer  established  a  1-gram  limit  to  the  out¬ 
put  charge  of  any  model  explosive.  Adherence  to  this  condition  was  part  of 
the  previously  mentioned  license  (Reference  6).  Accordingly,  the  test  condi¬ 
tions  under  which  each  of  the  bombs  could  be  simulated  are  the  particular 
combinations  of  accelerations  and  mass  lying  above  the  bold  line  in  Table  3. 
Table  4  lists  the  theoretical  dimensions  for  the  model  charges  listed  above 
the  bold  line  in  Table  3. 

A  review  of  the  results  presented  in  Tables  3  and  4  leads  to  the  conclu¬ 
sion  that  these  data  are  somewhat  representative  of  commercially  available 
detonators,  such  as  the  Reynolds  RP-83  detonator  shown  in  Figure  3.  This 
standard  detonator  consists  of  an  exploding  wire  bridge  (1),  a  low  density 
pressing  of  j«ntaerytRritoI  tetranitrate  (PETN)  (2),  a  high-density  RDX  ini¬ 
tiator  (3),  and  a  high-density  RDX  output  charge  (4),  all  of  which  are  con¬ 
tained  in  a  0.007-inch  thick  aluminum  cup  (5).  The  high-density  output  charge 


TABLE  4.  THEORETICAL  DIMENSIONS  OF  MODEL  EXPLOSIVES. 


Bomb 

Designation, 

pound 

60  g 

80  g 

100  g 

Diameter,  inch 

250 

0.17 

0.13 

0.10 

500 

0.23 

0.17 

0.14 

1000 

— 

0.22 

0.18 

2000 

— 

— 

0.22 

Height,  inch 

250 

0.45 

0.34 

0.27 

500 

0.54 

0.41 

0.33 

1000 

— — 

0.50 

0.40 

2000 

— 

1 

1 

1 

1 

0.51 

Case  Thickness, 

inch 

250 

0.005 

0.003 

0.003 

500 

0.005 

0.004 

0.003 

1000 

0.006 

0.005 

2000 

0.005 

Plastic  header 


Wire  bridge 


Aluminum  cup 
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Output  charge  RDX 


Figure  3.  Reynolds  Industries  RP-83  Detonator. 
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consists  of  four  individual  RDX  pressings  each  weighing  0.220  grams.  A  review 
of  these  features  indicated  that  if  the  length  of  the  aluminum  cup  could  be 
altered  by  the  manufacturer,  then  model  charges  having  the  following  charac¬ 
teristics  might  be  readily  obtained: 


Description 

Std. 

Mod.  1 

Mod.  2 

Mod.  3 

Mass,  gram 

0.880 

0.660 

0.440 

0.220 

Height,  inch 

0.700 

0.525 

0.350 

0.175 

Diameter,  inch 

0.275 

0.275 

0.275 

0.275 

Case  Thickness,  inch 

0.005 

0.005 

0.005 

0.005 

The  0.880-gram  detonator  is  the  standard  RP-83  detonator  and  it  required  only 
an  alteration  of  case  thickness  (Table  4).  Subsequent  discussions  with  the 
engineering  staff  of  Reynolds  Industries  led  to  a  contract  to  have  60  detona¬ 
tors  (15  of  each)  fabricated  for  use  on  this  project,  at  an  average  cost  of 
$11.63  each,  approximately  $2.20  above  the  cost  of  the  standard  RP-83.  Fig¬ 
ure  4  is  a  photograph  of  these  detonators.  An  additional  consideration  in  the 
selection  of  the  RP-83  was  the  safety  feature  of  the  exploding  wire  bridge 
(EW8).  The  EWB  is  virtually  insensitive  to  low  voltage  and  it  will  not  deto¬ 
nate  without  the  application  of  a  high  voltage;  thus  it  could  not  be  accident¬ 
ally  exploded. 

A  review  of  the  scaled  weights  for  the  weapons  listed  in  Table  3  and  the 
weights  of  the  model  charges  indicates  several  possibilities  for  exact  mass 
simulations  of  the  weapons  listed.  For  example,  the  yield  of  a  250-pound 
weapon  could  be  simulated  by 

0.220-gram  at  76.25  G 
0.440-gram  at  60.53  G 
0.660-gram  at  52.89  G 
0.880-gram  at  48.06  G 

Likewise,  the  yield  of  a  1000-pound  weapon  could  be  simulated  by 

0.440-gram  at  98.58  G 
0.660-gram  at  86.13  G 
0.880-gram  at  78.26  G 


However,  in  respect  to  characteristic  lengths  some  violation  of  similitude  had 
to  be  tolerated.  For  example,  a  review  of  Table  4  indicates  a  need  for  case 
thicknesses  ranging  from  0.003  to  0.006  inch.  Manufacturing  processes  limit 
the  case  thickness  to  a  minimum  value  of  0.005  inch.  This  standard  value  was 
selected  to  minimize  costs.  This  consideration  was  also  influenced  by  the 
need  to  have  similitude  in  the  intensive  properties  of  the  case  materials 
between  the  prototype  and  model.  The  former  values  were  essentially  unknown. 
Furthermore,  this  level  of  detail  was  considered  a  tertiary  consideration  in 
view  of  the  fact  that  both  the  prototype  and  model  configurations  stipulated 
that  the  weapon  would  be  embedded  in  the  concrete  burster  slab.  Accordingly, 
it  was  reasoned  that  the  confining  effect  of  the  burster  slab  would  influence 
the  free-field  blast  parameters  more  significantly  than  case  thickness. 

The  diameter  of  the  model  charges  was  a  uniform  0.250  inch.  This  offers 
a  more  serious  violation  of  similitude  as  indicated  in  Table  4  where  an  aver¬ 
age  value  of  0.16  inch  would  appear  more  appropriate.  The  effects  of  this 
demarcation  are  not  fully  known.  It  is  believed  that  the  effect  on  the 
recorded  peak  normal  stress  is  insignificant,  perhaps  affecting  attenuation  as 
reflected  in  the  duration  of  the  positive  phase  of  the  pressure  pulse.  Model 
charge  heights  were  also  standard,  being  multiples  of  0.175  inch.  Standard 
heights  and  diameter  were  necessary  for  economic  reasons;  alterations  to  the 
standard  0.250-inch  diameter  by  0.175-inch  high  pressing  were  cost  prohibi¬ 
tive.  Numerically,  the  violation  in  similitude  in  respect  to  height  is  less 
severe  than  that  for  diametrical  simulation.  The  effect  of  a  height  violation 
was  initially  judged  to  be  one  of  loading  the  burster  slab  to  alter  the  shape 
of  the  wave  front  in  the  sand  continuum.  Subsequent  analysis  of  the  data 
confirmed  this  judgement.  This  is  discussed  in  detail  in  Section  IV. 

The  model  charges  (Reynolds  detonators)  were  epoxied  in  place  in  the 
model  burster  slab  as  indicated  in  Figure  5.  The  detonators  were  epoxied  in 
place  to  simulate  the  condition  of  the  prototype  weapon;  that  is,  the  real 
weapon  does  not  fully  penetrate  the  burster  slab.  Under  this  condition  it  was 
assumed  that  the  real  weapon  would  firmly  lodge  in  the  burster  slab  prior  to 
detonation.  The  use  of  epoxy  to  bind  the  model  charge  (detonator)  to  the 
smooth  wall  of  the  model  burster  slab  was  an  attempt  to  simulate  the  real 
weapon  firmly  lodged  in  the  real  burster  slab.  The  epoxy  operation  was  per¬ 
formed  by  the  ordnance  technician  after  the  technician  took  charge  as  test 
conductor,  that  is,  after  the  test  fixture  containing  the  burster  slab  and 
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sand  was  attached  to  the  swing-bucket  platform  on  the  centrifuge.  After  the 
epoxy  was  set,  and  with  the  centrifuge  rotating  at  the  prescribed  angular 
velocity,  the  charge  was  detonated.  A  Reynolds  Industries  Model  FS-10  control 
panel  and  firing  module  were  used  to  initiate  the  EWB  within  the  detonator 
(model  charge).  The  control  panel  was  mounted  externally  to  the  centrifuge, 
while  the  firing  module  was  mounted  atop  the  vertical  axis  of  the  centrifuge. 
Voltage  was  applied  to  the  firing  module  via  a  coaxial  cable  and  a  slipring, 
both  of  which  were  physically  isolated  from  the  data  acquisition  cables  and 
sliprings.  The  capacitor  discharge  firing  module  produced  3000  volts  d.c. 
However,  even  with  the  physical  separation  of  this  system  from  the  data 
acquisition  system,  both  systems  were  not  fully  isolated.  Thus  the  high 
frequency  associated  with  the  capacity  discharge  of  the  FS-10  was  recorded  on 
the  data  channels. 


EXPERIMENTAL  PROGRAM 

Equation  (1)  is  rewritten  here  for  ease  of  discussion. 
,1/3 


,1/3 


-o J  pc2  C2  C2  a2  c2  i  P0 

If  Equation  (1)  is  considered  to  model  the  blast-induced  normal  stress, 
a,  in  the  centrifuge,  that  is,  if  material  properties  do  not  vary  from  test  to 
test  and  if  characteristic  lengths  are  scaled  according  to  the  G-level,  then 
for  the  centrifuge  environment.  Equation  (1)  becomes 
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If  the  same  type  of  explosive  is  used  for  all  tests,  then  the  dependence 
on  P0  is  eliminated  and  the  energy  of  the  explosive  can  then  be  replaced  with 
its  mass.  Thus  a  simpler  form  of  Equation  (1)  becomes 
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where  the  characteristic  length  is  replaced  with  R,  the  range  at  which  a  is  to 
be  experimentally  evaluated.  The  terms  that  remain  in  Equation  (1)  are  not 
nondimensional ;  that  is,  they  are  not  true  "-terms.  This  does  not  invalidate 
Equation  (2);  rather,  it  simply  means  that  the  magnitude  of  each  term  has  been 


altered;  each  remains  fundamental  to  the  physical  process  being  modeled.  If 
it  is  assumed  that  Equation  (1)  is  complete,  it  can  be  written  as 


Equation  (3)  defines  a  surface  in  a  three-dimensional  system  as  shown  in 
Figure  6,  where  the  G-level  is  confined  to  the  energy  term,  G(W)1/3.  Experi¬ 
mental  data  for  the  scaled  distance  R(W_1/3)  should  not  reveal  any  gravity 
affects  because  distance  scales  as  1/n  and  mass  as  1/n3.  Mathematically,  at 
least,  the  scaling  factors  are  thus  equal  for  both  the  numerator  and  denomina¬ 
tor,  and  cancel  each  other.  The  purpose  of  the  experiments  was  to  confirm 
this  observation  and  to  demonstrate  that  centrifugal  blast-induced  stress  data 
do  scale  as  indicated  by  Equation  (3). 

To  demonstrate  scalability,  12  experiments  were  planned  as  outlined  in 
Table  5.  In  reviewing  Table  5, it  is  important  to  note  that  the  normal  stress 
measurements  were  taken  at  the  same  scaled  range,  but  that  the  gages  were  not 
at  the  same  range.  For  example,  consider  the  series  of  four  tests  conducted 
at  50  G  and  X  =  2: 


V  m9 

Wp,  lb 

R,  in 

220 

3.06 

1.47 

440 

3.86 

1.86 

660 

4.42 

2.12 

880 

4.86 

2.33 

These  data  demonstrate  the  term  scalability:  adjustment  of  the  G-level  to 
scale  the  yield  to  produce  scaled  distances  having  the  same  value,  but  derived 
from  different  explosive  weights  and  radial  offsets,  R.  Radial  offsets, 
standoff  distances,  were  located  with  respect  to  the  bottom  of  the  explosive 
charge  (see  Figure  5). 

MOOEL  CONSIDERATIONS 

Modeling  of  the  Burster  Slab 

For  purposes  of  modeling,  the  prototype  burster  slab  was  considered  to  be 
48  inches  thick  with  No.  4  reinforcing  bars  (temperature)  at  14  inches,  center 
to  center,  both  ways,  in  both  the  top  and  bottom  of  the  slab.  Rather  than 
model  the  slab  thickness  exactly  for  each  G-level  used  in  the  experiments, 
model  slab  thicknesses  were  established  as  indicated  in  Table  6. 


<0 


(U 


CO 

t- 

3 

JD 


U) 

<u 


s- 

0) 

CO 


T3 

o> 

+-> 

u 


o 

o 


o  • 

i/>  -*-> 
• —  CO 
JO  0) 


to 

JO  0) 

2  J= 


CO 

<u 


L. 

Q)  "O 
co  o> 
“O 
CO  3 


x:  u 

h-  C 
JO  -r- 


tHTT 


TABLE  6.  MODEL  BURSTER  SLAB  THICKNESSES. 


— 

G 

Simul ation 
Thickness, 
inches 

Selected 

Thickness, 

inches 

31.44 

1.53 

1.40 

39.60 

1.21 

1.40 

46.37 

1.04 

0.97 

50.00 

0.96 

0.97 

53.00 

0.91 

0.97 

71.88 

0.67 

0.64 

79.52 

0.60 

0.64 

The  temperature-reinforcing  steel  was  modeled  using  standard  grade  hard¬ 
ware  cloth  made  of  23-gage  wire.  The  hardware  cloth  has  four  meshes  per  lin¬ 
ear  inch  and  is  commercially  available,  but  only  on  special  order.  In  respect 
to  simulation,  it  accurately  models  the  prototype  steel  percentage  at  a  scale 
factor  corresponding  to  50  G,  the  median  of  the  G-levels  reached  in  these 
experiments. 

The  concrete  was  simulated  using  a  laboratory-developed  grout.  Labora¬ 
tory  trial  mixes  included  determination  of  unit  weights  and  7-day  compressive 
strengths.  The  selected  mix  had  the  properties  listed  in  Table  7. 

Circular  burster-sl abs  21  inches  in  diameter  were  cast  from  mixes  having 
the  proportions  stated  in  Table  7.  Test  cylinders  were  also  cast  from  each 
mix.  Each  model  slab  was  encased  in  plastic  for  a  cure  period  of  3  days 
before  removal  from  its  circular  steel  form.  It  was  then  placed  in  a  moist 
room  and  cured  for  at  least  28  days  before  being  used  in  a  centrifugal  test. 
Figure  7  is  a  postshot  view  of  a  model  burster  slab.  The  wires  extending  from 
the  slab  on  each  side  of  the  crater  were  cast  into  each  slab  to  facilitate  the 
handling  of  the  slab. 
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TABLE  7.  PROPERTIES  OF  THE  MODEL  BURSTER  SLAB 
SIMULATED  CONCRETE  MIX. 


Cement,  Type  I 

1.0 

Water-Cement  Ratio 

0.6 

Sand 

2.33 

Baroid 

1.0 

Plasticizer 

0.03 

Unit  Weight 

147  1 b/f t 3 

Compressive  Strength,  28-days 

3300  lb/ in2 

As  indicated  in  the  introduction  to  this  section,  the  burster  slab  is 
supported  by  a  sand  continuun.  No  specifications  were  provided  detailing  the 
engineering  properties  of  the  prototype  sand,  other  than  to  assume  it  to  be  a 
quality  sand  such  as  that  used  in  a  portland  cement  concrete  mix.  Such  a 
material  would  have  a  gradation  ranging  from  the  No.  4  to  No.  200  sieve  as 
indicated  in  Figure  8.  Application  of  the  scaling  relationships  suggests  that 
a  sand  used  in  the  model  should  have  the  particle  size  distribution  indicated 
on  the  extreme  right  of  Figure  8.  At  the  G-levels  attained  in  this  study,  a 
strict  application  of  the  scaling  relationship  indicates  a  need  for  a  material 
with  a  particle  size  distribution  ranging  from  the  No.  200  size  to  the  1- 
micron  size.  Particles  within  that  size  range  begin  to  display  clay-like 
behavior,  a  significant  departure  from  the  granular  behavior  of  the  prototype 
system.  Additionally,  strict  application  of  the  concepts  of  similitude  would 
require  a  separate  gradation  for  each  G-level  and,  more  importantly,  consider¬ 
ation  of  strain-rate  affects.  This  latter  consideration  is  not  well  under¬ 
stood  and  the  former  suggests  the  expenditure  of  considerable  effort  to  pre¬ 
pare  each  gradation.  Because  of  these  considerations  it  was  decided  to  use  a 
commercially  available,  fine-grained,  well-graded  sand;  its  gradation  is  also 
shown  on  Figure  8.  This  sand  scales  the  prototype  material  in  the  10-  to  15-G 
range.  It  was  used  in  all  the  tests,  regardless  of  the  G-level. 

Crystal  Silica  No.  70  sand  has  the  following  engineering  properties: 

Specific  Gravity  2.65 

Relative  Density  ASTM  D-2049 

Minimum  74  lb/ft3 

Maximum  94  lb/ft3 

Angularity  (Reference  7)  Angular 

A  primary  consideration  for  support  of  the  test  program  was  related  to 
the  consistent  production  of  uniform  sand  models.  The  customary  technique  of 
raining  sand  into  the  test  fixture  was  attempted,  but  this  procedure  failed  to 
produce  models  having  a  consistently  uniform  density.  The  raining  technique 
apparently  failed  because  the  sand  was  very  fine  grain.  Several  other  tech¬ 
niques  were  attempted  without  satisfying  the  placement  criterion.  A  vibration 
procedure  was  finally  adopted. 


Figure  8.  Grain  Size  Distribution  Graph. 


Two  Airoviber  vibrators  were  rigidly  attached  to  the  base  of  the  test 
fixture  (Figure  9).  These  air-driven  vibrators  produce  800  pounds  of  centrif¬ 
ugal  force  each,  when  subject  to  compressed  air  at  80-lb/in2  gage.  The  sand 
was  loosely  placed  in  the  fixture  in  2-inch  layers  and  vibrated  for  10  min¬ 
utes.  This  procedure  was  repeated  until  the  sand  was  10  inches  deep  in  the 
fixture.  Actual  sand  densities  were  based  upon  gross  weight  and  the  volume 
occupied  by  the  sand.  The  mean  density  for  all  tests  was  103.9  lb/ft3,  the 
range  101.9  to  105.9  lb/ft3,  the  standard  deviation  0.80  lb/ft3,  and  the 
coefficient  of  variation  7  percent. 

DATA  ACQUISITION 

Stress  Gage  Development 

Normal  stresses  within  the  sand,  generated  by  the  explosion  of  the  model 
charge,  were  measured  using  a  foil  pressure  sensor.  Dynasen's  carbon-foil 
pressure  sensors,  C300-50-EKRTE,  were  used  in  this  study.*  These  sensors  were 
selected  primarily  for  size  considerations,  having  nominal  dimensions  of  0.050 
by  0.060  by  0.0025  inch.  At  the  maximum  G-level  of  the  centrifuge  these 
dimensions  scale  to  5  by  6  inches,  small  enough  to  represent  a  point  measure¬ 
ment  in  the  prototype  field.  However,  all  of  the  experiments  were  performed 
at  G-levels  below  this  value;  thus  the  prototype  measurement  area  was  smaller, 
perhaps  on  the  order  of  4  square  inches  on  the  average.  An  exploded  and 
enlarged  view  of  the  Dynasen  pressure  sensor  is  shown  as  Figure  10. 

By  necessity  these  sensors  had  to  be  attached  to  a  mounting  medium  in 
order  to  place  the  sensor  in  the  sand  continuum.  Having  no  prior  experience 
with  these  sensors  in  a  free-field  application,  several  sensor  mounting  tech¬ 
niques  were  tried  before  a  successful  gage  was  developed.  Thus  the  initial 
tests  served  as  a  gage-development  program.  Only  several  acceptable  data 
records  were  obtained  from  these  tests  because  o^  ^--h  problems  as  debonding, 
frequency  response,  debris  impact,  and  strain  effects.  The  gage  that  eventu¬ 
ally  evolved  consisted  of  the  carbon-foil  pressure  sensor  bonded  to  0.008-inch 
stainless  shim  stock,  cut  from  sheet  rather  than  rolled  shim  stock.  The  shim 
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stock  measured  nominally  2.5  inches  long  by  0.5  inch  wide.  The  length  was 
necessary  to  accommodate  the  2-inch  foil  tabs  which  were  attached  to  the 
carbon-foil  sensor.  The  sensor  and  its  tabs  were  bonded  to  the  shim  stock 
with  EPY-150.  The  tabs  of  the  pressure  sensor  were  connected  to  a  4-wire  1122 
Alpha,  forming  the  pressure  gage  shown  in  Figure  11.  These  gages  were  pre¬ 
cisely  placed  in  the  sand  continuum  as  suggested  in  Figure  11. 

Instrumentation  System 

A  major  problem  with  the  centrifuge  cabling  system,  as  it  existed  prior 
to  this  study,  was  that  the  28  data  channels  were  made  from  single-conductor 
coaxial  cables  with  unshielded  terminals.  Pretest  evaluation  indicated  the 
performance  of  the  system  to  be  poor.  It  was  characterized  by  excessive 
cross-talk  and  a  variable  frequency  response  (depending  on  the  termination 
resistance)  with  a  maximum  frequency  response  of  10  MHz.  To  improve  the 
response  of  this  system,  six  of  the  data  cables  were  replaced  with  Belden- 
type,  18-gage,  twisted-pair  cable,  each  having  a  maximum  length  of  40  feet. 
Based  upon  experience  with  systems  of  similar  characteristics,  the  frequency 
response  of  the  upgraded  system  was  estimated  to  exceed  the  maximum  frequency 
of  interest  (240  MHz).  However,  it  was  impossible  to  shield  the  slip-rings  as 
well  as  the  data  cables  on  the  boom  of  the  centrifuge.  These  deficiencies  led 
to  fire-set  noise  on  the  data  traces  and  made  some  records  difficult  to 
interpret. 

Figure  12  is  a  wiring  diagram  for  the  stress  gages.  Figure  13  is  a  block 
flow  diagram  of  the  instrumentation  system.  All  data,  including  calibration 
and  fiducial,  were  recorded  on  standard  wide  band  FM440.  Playback  was  on  a 
Nicolet  digital  oscilloscope  with  a  hard  copy  from  an  x-y  recorder. 

STRESS  GAGE  PERFORMANCE 

Eighty-two  pressure  time  histories  were  obtained  from  the  series  of  tests 
which  included  a  burster  slab.  These  records  were  examined  for  determination 
of  peak  normal  stress,  rise  time,  stress  duration,  and  time  of  arrival.  This 
included  records  obtained  during  the  early  phase  of  the  testing  program  which 
was  essentially  a  gage-development  phase.  All  of  the  data  were  not  considered 
valid  because  of  mechanical  failure  of  the  gage,  noise  on  the  record,  or  mal¬ 
function  of  the  explosive.  Thirty  of  the  82  data  points  were  considered  to  be 
invalid.  Figure  14  is  a  pictorial  representation  of  the  observed  responses  of 


Figure  12.  Configuration  Detail— Dynasen  Shock  Pressure  Gage. 


the  gage  which  were  considered  to  be  unacceptable.  The  mineral  in  the  upper 
corner  of  each  frame  reflects  the  number  of  data  traces  rejected  for  the  rea¬ 
son  indicated.  Thirteen  records  from  the  first  three  tests  (all  repeats  of 
the  660-mg  at  50  G)  were  rejected,  principally  because  of  strain  effects  and 
broken  gages.  The  other  17  rejected  records  were  scattered  uniformly  through¬ 
out  the  remaining  11  tests.  Seven  of  these  17  records  were  rejected  because 
of  debonding,  the  negative/positive  rejection  mode. 

Twenty-three  pressure  gages  were  used  in  the  test  series  conducted  with¬ 
out  a  burster  slab;  all  of  these  records  were  considered  acceptable. 

Figures  15  and  16  are  selected  normal  stress/time  records  obtained  in 
this  study. 


SECTION  IV 
RESULTS 


TESTS  WITH  A  BURSTER  SLAB 

Table  8  is  a  summary  of  the  data  obtained  from  those  tests  in  which  a 
burster  slab  was  included  in  the  test  configuration.  These  data  include  52 
recordings  of  normal  stress,  a,  taken  at  scaled  distances,  X,  ranging  from  0.5 
to  5.  The  units  are  lb/in2  and  ft/W1/3  respectively.  The  data  also  include 
the  sand  dilatational  wave  velocity,  C2  (ft/s).  This  is  a  calculated  result 
based  upon  the  known  gage  location,  R,  and  a  determination  of  the  time  of 
arrival  of  the  stress  wave  at  each  gage.  These  data  were  plotted  on  normal 
probability  graph  paper  to  determine  whether  the  distribution  follows  the 
pattern  of  a  normal  curve.  These  plots  are  presented  as  Figures  17,  18,  and 
19.  Rather  than  plotting  the  dependent  variable,  cr,  these  data  are  presented 
in  the  nondimensional  form  o/p2C \  which  is  considered  to  be  a  more  fundamen¬ 
tal  term  for  ground  motion  studies  than  a. 

A  review  of  these  figures  indicates  a  tendency  for  the  data  to  plot  in  a 
linear  manner.  This  linearity,  though  subjective,  suggests  that  the  data 
follow  the  pattern  of  a  normal  curve.  The  implication  is  that  the  data  con¬ 
tain  primarily  only  random  error.  In  an  attempt  to  control  experimental 
errors,  all  test  models  (placement  of  sand  and  pressure  gages)  were  prepared 
by  the  same  technician.  In  an  effort  to  minimize  the  effect  of  improved 
skills  (systematic  error),  each  test  (a  different  combination  of  model  charge 
weight  and  G-level)  was  assigned  a  test  number  derived  from  a  random  draw. 

A  review  of  the  normal  probability  plots  indicates  an  average  coefficient 
of  variation  of  13  percent  for  all  x-groups.  The  range  in  this  statistic  is  9 
to  20  percent.  When  examined  in  terms  of  normal  stress,  the  data  indicate  an 
average  variation  about  the  mean  of  7  percent  at  one  standard  deviation.  The 
average  range  in  normal  stress  about  the  mean  value  is  23  percent.  These 
statistics  indicate  that  the  quality  control  for  the  experiments  was  excellent 
for  this  type  of  study,  and  that  the  experiments  contained  a  high  degree  of 
replication. 


Normal  Probability  Plot  of  Nondimensional  Stress 
for  Tests  with  a  Burster  Slab,  X  =  0.5  and  0.8. 


Figure  19.  Normal  Probability  Plot  of  Nondimensional  Stress 
for  Test  with  a  Burster  Slab,  X  =  3  and  5. 


All  sources  of  error  inherent  to  the  experiments  are  not  known.  A  gage 
or  explosive  charge  placement  error  of  0.02  inch  (a  reasonable  expectation) 
results  in  an  average  error  of  4  percent  in  the  scaled  distance.  The  close-in 
gages  are  more  adversely  affected  by  placement  errors  than  those  gages  at  a 
larger  range.  However,  no  statistic  could  be  found  that  suggests  other  than  a 
random  error  with  respect  to  gage  or  explosive  location.  Pressure  sensor 
calibration  was  provided  by  the  manufacturer.  Each  sensor  could  not  be  cali¬ 
brated  because  the  process  is  essentially  destructive  to  the  sensor.  The 
sensors  are  manufactured  in  small  batches.  Sensors  randomly  selected  from 
each  batch  are  calibrated.  The  manufacturer  stipulates  that  the  maximum  vari¬ 
ation  for  sensors  from  different  batches  is  5  percent  from  the  calibration 
provided  and  less  than  5  percent  intrabatch.  However,  the  sensors  were 
mounted  on  shim  stock  to  form  a  stress  gage.  The  Civil  Engineering  Research 
Facility  (CERF)  did  not,  however,  have  the  equipment  required  to  calibrate 
these  gages,  so  the  manufacturer's  calibration  had  to  be  used. 

Probably  the  largest  single  source  of  random  error  lies  in  the  determina¬ 
tion  of  time  of  arrival  used  to  compute  dilatational  velocity,  C2.  An  error 
in  the  determination  of  the  time  of  arrival  on  the  order  of  1  microsecond  (a 
reasonable  expectation)  results  in  a  variation  of  the  dilatational  velocity  of 
nearly  4  percent.  This  translates  into  an  error  in  p/p2C2  of  more  than 
8  percent.  The  accumulation  of  all  of  these  known  sources  of  error  easily 
exceeds  the  variation  observed  in  the  data.  Clearly,  the  quality  control  on 
the  experiment  was  acceptable. 

The  variable  of  primary  interest  in  this  study  is  the  normal  stress, 
rather  than  the  ratio  c/p2C|.  The  ultimate  application  is  that  of  pre¬ 
dicting  o  as  a  function  of  the  scaled  distance,  X.  Therefore  the  paired 
observations  of  these  variables  were  subjected  to  multiple  linear  regression 
analyses  to  determine  the  functional  form  of  the  best-fit  relationship  between 
o  and  X.  The  criterion  of  least  squares  was  used  to  provide  an  estimate  of 
the  equation  of  the  regression  line.  The  strength  of  the  fit  was  evaluated  on 
the  basis  of  the  correlation  coefficient  calculated  for  each  linear  relation¬ 
ship  obtained  from  the  analysis.  This  statistic  was  selected  because  the  data 
points  (o^,  X.j )  for  i  =  1,  2.. 52,  are  values  of  a  pair  of  random  vari¬ 
ables.  That  is,  neither  variable  can  be  assumed  to  be  known  without  error. 


and  thus  the  data  are  a  sample  of  a  bivariate  population.  The  assunption  that 
the  bivariate  sample  has  a  normal  distribution  (the  errors  are  random)  is 
derived  from  Figures  17,  18,  and  19,  and  the  foregoing  discussion  concerned 
with  experimental  errors. 

Figure  20  is  a  plot  of  the  paired  values  of  a^,  X^  presented  in 
Table  8.  These  data  are  plotted  on  log-log  scales.  This  transformation  per¬ 
mits  a  linear  regression  analysis  and  was  necessary  because  a  suitable  fit  to 
the  data  could  not  be  found  in  the  arithmetic  domain.  A  linear  relationship 
is  necessary  for  a  least-square  fit.  Least  squares  estimators  of  the  true 
relationship  have  the  smallest  variance  and  are  the  most  reliable  in  the  sense 
that  they  are  subject  to  the  smallest  chance  variation.  The  best  single¬ 
variable  model  fit  to  the  data  is 

a  *  5879X-1*325  (4) 

where 

o  =  vertical  normal  stress,  lb/in2 
x  =  scaled  distance,  ft/1 b1 7 3 

for  vihich  the  sample  correlation  coefficient,  r2,  has  a  value  of  0.979.  This 
value  indicates  that  98.9  percent  of  the  sample  variation  in  log  a  can  be 
attributed  to  the  linear  variation  in  log  X.  By  application  of  the  Fisher 
transformation,  the  95-percent  confidence  interval  of  0.982  <  r  <  0.994  has 
been  established  as  an  estimate  of  the  true  strength  of  the  linear  relation¬ 
ship  between  log  a  and  log  X.  The  standard  error  of  estimate  in  the  exponent 
is  0.027;  for  the  coefficient,  it  is  0.01017  applied  to  its  logarithm.  Equa¬ 
tion  4  is  plotted  on  Figure  20.  Evidence  of  the  goodness  of  fit  is  revealed 
in  Figure  21.  This  figure  is  a  plot  of  the  residual  (log  ck  -log  a^)  versus 
log  X.. ,  where  oi  is  the  predicted  value  of  The  inference  that  the  fit  is 
good  is  based  upon  the  observation  that  the  residuals  appear  to  be  equally 
distributed  about  the  zero  residual. 

Other  multiple-variable,  linear  models  were  developed  in  an  attempt  to 
improve  on  the  correlation  coefficient.  However,  no  significant  improvement 
in  this  statistic  was  detected.  A  summary  of  these  models  follows. 


log  o  =  -1.325  log  X  +  5.9277 


r2  =  0.9792 


a  =  normal  stress,  lb/in2 
X  =  scaled  distance,  ft/1b1/3 

log  a  =  -1.325  log  X  -  0.098  log  X'  +  5.9670  r2  =  0.9794 

X'  =  scaled  depth  of  burster  slab 

log  o  =  -1.327  log  X  -  0.030  log  X"  +  5.9279  r2  =  0.9794 

X"  =  scaled  charge  length 

log  o  =  -1.324  log  X  -  0.101  log  X'  -  0.031  log  X"  +  5.9715  r2  =  0.9796 

In  reviewing  Figure  20  it  is  important  to  recognize  that  the  data  points 
at  a  particular  scaled  distance  contain  only  a  few  repeat  tests.  To  the  con¬ 
trary,  the  ordinate  values  are  derived  from  tests  with  explosive  yields  scaled 
according  to  the  -simi 1 itude  relationship  developed  in  Section  II.  The  range 
in  data  is  due  to  experimental  error.  None  of  the  forms  of  the  dependent 
variable  (o1 ,  a/P0,  a/p2C2)  revealed  significant  correlation  with  the  G- 
level;  in  a  regression  analysis,  the  maximum  observed  value  of  the  correlation 
coefficient  was  0.144.  Thus  the  data  in  the  o,X-plane  (Figure  20)  does  not 
contain  any  gravity  effects.  These  data  were  also  analyzed  to  determine  if 
any  consistent  trends  could  be  observed  relative  to  the  weight  of  the  explo¬ 
sive,  V^.  The  data  were  grouped  according  to  the  weight  of  the  explosive; 
and  then  a  linear  regression  was  performed  to  determine  the  coefficients  in 
the  linear  relationship  between  o/p2C2  and  X.  This  form  of  the  dependent 
variable  was  chosen  because  it  contains  a  larger  random  error  than  a.  The 
results,  shown  in  Table  9,  indicate  no  significant  difference  between  the 
individual  groups  and  that  for  the  entire  group  of  weights.  There  is  an  indi¬ 
cation  that  the  attenuation  decreases  with  increasing  weight.  The  effect  is 
not  due  to  weight;  it  is  due  to  shape  of  the  charge.  The  shape  of  the  model 
charge  becomes  increasingly  cylindrical  as  its  weight  increases  (see 
Section  III). 


E' 


w 

"m 

No.  of  Tests 

Equation3 

r2 

220  mg 

18 

o/p2C2  =  0.389  X-1'219 

0.980 

440  mg 

5 

o/ p2C2  =  0.318  X-l-386 

0.997 

660  mg 

9 

a/p2C2  =  0.359  X-1  -304 

0.985 

880  mg 

20 

o/p2C2  =  0.360  X-1 *292 
*•  2 

0.958 

Combined 

52 

a/p2C2  =  0.357  X-1*271 
*  2 

0.970 

Feet-slugs-seconds 

The  other  ir-term  involved  in  these  experiments  is  G(W)1/3.  Because  W  in 
that  term  is  (i^  •  G3) ,  no  gravity  effects  should  be  detected  in  a  linear 
regression  analysis  with  o  as  the  dependent  variable  and  G(W)1/3  as  the  inde¬ 
pendent  variable,  at  any  value  of  X.  Thus,  because  G(W)1/3  is  an  energy  term, 
at  X  =  0  the  detonation  pressure  is  independent  of  the  energy  and  the  o  versus 
G(W)1/3  relationship  should  have  a  zero  slope.  Likewise,  at  any  other  value 
of  X,  the  slope  should  remain  zero,  except  that  the  value  of  o  decreases  with 
increasing  x.  The  relationship  is  the  surface  previously  shown  in  Figure  6. 

Linear  regression  analyses  of  the  data  in  Table  8  confirmed  the  shape  of 
the  surface.  A  linear  function  of  the  form  a  =  m[G(W)1/3]  +  b,  was  fit  to  the 
data  at  X  =  0.5,  1,  2,  and  5.  In  all  cases,  except  X  *  0.5,  the  slope  values 
(m)  were  such  that  -0.4  <  m  <  4  and  the  intercept  values  (b)  were  in  accept¬ 
able  agreement  with  the  stress  predicted  by  Equation(4)at  the  corresponding 
values  of  X. 


As  noted,  the  results  at  X  =  0.5  did  not  agree  with  the  assumed  shape  of 
Figure  6.  The  best-fit  line  for  X  *  0.5  had  the  form  a  =  -17[G(W)1/3]  + 
21,394,  with  a  correlation  coefficient  of  0.44.  While  this  indicates  a  low 
correlation,  the  slope  departs  significantly  from  the  former  values  (-0.4  < 
m  <  4)  and  the  intercept  is  nearly  40  percent  larger  than  the  value  given  by 
Equation(4)at  X  *  0.5.  A  review  of  the  data  for  X  =  0.5  indicates  that  the 


departure  from  a  zero  slope  is  due  to  an  error  in  the  design  of  the  experi¬ 
ment.  As  it  developed,  the  heavier  charges  (880  mg)  were  accelerated  to  the 
lowest  G-levels,  while  the  smallest  (220  mg)  were  taken  to  high  G-levels.  The 
result  is  that  at  low  values  of  G(W)1/3  the  charges  were  truly  cylindrical, 
while  at  high  values  of  G(W)1/3  the  charge  shape  was  more  nearly  spherical. 
Hence,  a  case  of  cylindrical  divergence  versus  spherical  divergence,  or  higher 
stresses  close  in  for  the  880  mg  charges.  A  tabulation  of  selected  data 
clearly  shows  the  trend. 


G(W) 1/3 

Wm 

m 

17,700  lb/in2 

306 

880  mg 

13,730  lb/in2 

400 

660  mg 

12,000  lb/in2 

482 

220  mg 

Interestingly,  if  the  data  points  associated  with  the  880  mg  are  removed 
from  the  data  set  at  A  =  0.5,  the  slope  falls  into  the  range  reported  above 
(-0.4  <  m  <  4)  and  the  intercept  is  within  10  percent  of  that  predicted  by 
Equation(4).  If  all  the  data  are  fit  to  a  model,  based  upon  the  considerations 
of  cylindrical  divergence,  both  the  slope  and  intercept  remain  unacceptable, 
because  the  smaller  charges  do  not  simulate  cylindrical  shapes.  The  choice 
then  is  to  remove  the  220  mg  data  from  the  set;  the  fit  improves,  m  <  7,  and 
the  intercept  exceeds  the  stress  given  by  Equation(4)by  11  percent.  The  clear 
conclusion  is  that  at  A  <  0.5  charge  shape  is  an  important  factor. 

TESTS  WITHOUT  A  BURSTER-SLAB 

A  series  of  four  tests  without  a  burster  slab  was  performed  for  the  pri¬ 
mary  purpose  of  investigating  the  influence  of  the  burster  slab  on  the  attenu¬ 
ation  of  normal  vertical  stress,  o.  This  series  was  planned  after  the  results 
from  the  tests  having  a  burster  slab  indicated  less  attenuation  than  pre¬ 
dicted.  All  details  in  this  test  series  were  identical  with  those  tests  which 
included  a  burster  slab.  The  model  charge  was  placed  at  a  depth  in  the  sand 
corresponding  to  the  thickness  of  the  burster  slab,  had  it  been  present. 


Because  all  tests  in  this  series  were  performed  at  50  6  the  depth  of  burial 
was  a  constant  0.97  inch.  The  data  obtained  from  this  series  are  presented  in 
Table  10.  Radial  offsets,  standoff  distances,  were  located  with  respect  to 
the  bottom  of  the  explosive  charge.  Interestingly,  23  gages  were  used  in  this 
series  and  none  of  the  data  was  rejected  for  any  of  the  reasons  outlined  in 
Section  III.  The  normal  distribution  tendency  of  these  data  is  apparent  from 
a  review  of  Figures  22  and  23.  Figure  24  is  a  plot  of  these  data,  vertical 
stress  o  (lb/in2)  versus  scaled  distance  X  (ft/lb1/3). 


TABLE  10.  SUMMARY  OF  TEST  DATA  FOR  TESTS 
WITHOUT  BURSTER  SLAB. 


Observation 

X, 

ft/1 b1 /3 

w, 

mg 

G 

y2« 

lb/f tJ 

o, 

lb/in2 

1 

0.5 

220 

50 

103.44 

20070 

2 

0.5 

220 

50 

103.44 

16880 

3 

0.5 

440 

50 

105.40 

13610 

4 

0.5 

440 

50 

105.40 

11820 

5 

0.5 

660 

50 

103.54 

17130 

6 

0.5 

660 

50 

103.54 

16500 

7 

0.5 

880 

50 

103.60 

14200 

8 

1.0 

220 

50 

103.44 

1690 

9 

1.0 

220 

50 

103.44 

1950 

10 

1.0 

440 

50 

105.40 

1730 

11 

1.0 

440 

50 

105.40 

2070 

12 

1.0 

660 

50 

103.54 

1930 

13 

1.0 

660 

50 

103.54 

2110 

14 

1.0 

880 

50 

103.60 

1840 

15 

1.0 

880 

50 

103.60 

1810 

16 

2.0 

220 

50 

103.44 

375 

17 

2.0 

220 

50 

103.44 

240 

18 

2.0 

440 

50 

105.40 

390 

19 

2.0 

440 

50 

105.40 

453 

20 

2.0 

660 

50 

103.54 

511 

21 

2.0 

660 

50 

103.54 

630 

22 

2.0 

880 

50 

103.60 

305 

23 

2.0 

880 

50 

103.60 

235 

Figure  22.  Normality  Plot  of  Normal  Stress  for  Tests  Without 
Burster  Slab.  X  *  0.5  and  1. 


The  best-fit  equation  for  these  data  is 

a  =  2,100  X-2-684  (5) 

for  which  the  correlation  coefficient  is  0.99.  The  standard  error  of  estimate 
for  the  attenuation  factor  is  0.096  and  0.0234  for  the  logarithm  of  the 
intercept. 

A  discussion  of  the  results  presented  in  this  section  is  included  in 
Section  V. 


(The  reverse  of  this  page  is  blank.) 


SECTION  V 

DISCUSSION  OF  RESULTS 


FREE-FIELD,  NORMAL  STRESS 

A  plot  of  the  best-fit  line  for  the  normal  stress  scaled  distance  (a 
versus  X)  relationship,  for  the  tests  without  a  burster  slab,  is  shown  as 
Figure  25a.  The  line  labeled  Data  is  Equation  (5)  presented  in  the  former 

section.  _ this  line  is  the  predicted  preshot  relationship  which  was 

used  to  set  instrumentation  levels.  The  predicted  relationship  is  an  average 
of  predictions  derived  from  (1)  the  nonnuclear  design  manual  (p2C2  =  5,168  lb¬ 
s/ft3,  coupling  factor  =  0.7  TNT),  (2)  a  compilation  of  high-explosive  test 
events,  and  (3)  data  in  Reference  1.  Figure  25a  shows  good  agreement  between 
the  predictions  and  the  actual  measurements.*  Figure  25a  is  a  replot  of 
Figure  24.  A  review  of  Figure  24  reveals  that  the  best-fit  line  [Equa¬ 
tion  (5)]  passes  slightly  above  the  data  points  at  X  =  1.  This  suggests  that 
a  curvilinear  relationship  might  better  fit  the  data.  It  might  also  suggest 
that  a  break  (change  in  slope)  occurs  in  the  linear  relationship  so  that  the 
function  remains  linear,  but  becomes  steeper  at  x  <  0.5.  Such  a  trend  is 
well  known.  However,  concern  must  be  expressed  over  the  range  in  the  data  at 
x  *  2.  At  the  low  stress  levels  associated  with  X  =  2  the  pressure  sensors 
were  exposed  to  a  stress  on  the  order  of  0.05  percent  of  their  range  (50-kilo- 
bars),  and  nearly  0.5  kilobar  below  the  nearest  calibration  pressure.  How¬ 
ever,  real  concern  is  not  expressed  over  this  point,  because  the  tests  without 
a  burster  slab  were  performed  only  to  confirm  the  attenuation  determined  for 
the  tests  having  a  burster  slab  in  the  test  configuration. 

Figure  25b  is  a  plot  which  presents  a  comparison  between  the  results  for 
the  tests  without  a  burster  slab  (a  =  2100  x"2‘684)  and  the  results  for  the 
tests  with  a  burster  slab  (o  =  5879  x"1'325).  The  two  lines  shown  are  graph¬ 
ical  representations  of  these  best-fit  equations.  As  noted,  the  effect  of  the 


♦Items  (1)  and  (2)  were  provided  by  Applied  Research  Associates,  Incorporated 
of  Albuquerque,  New  Mexico.  This  firm  provided  general  consultation  to  the 
principal  investigator  on  this  study.  The  prediction  equation  is  only  valid 
for  lambda  between  2  and  15.  However,  because  no  known  data  existed  in  the 
range  of  0.5  <  X  <  2,  the  prediction  equation  was  arbitrarily  extended  into 
this  range  to  provide  an  estimate  of  anticipated  normal  stress  for  stress  gage 
calibration,  center  frequencies,  and  band-edge  calibrations. 


Figure  25.  Comparison  of  Test  Results,  Predictions,  and  SAMSON  Code  Values 


burster  slab  is  that  of  increased  stress  at  X  >  0.5.  Apparently,  the  higher 
seismic  wave  speed  in  the  burster  slab  creates  more  of  a  planar  shock  front  in 
the  sand  as  suggested  by  a  -1.325  attenuation  factor. 


A  series  of  two-dimensional  calculations  was  performed  to  support  the 
conclusion  drawn  from  the  tests  with  a  burster  slab  present;  namely,  the 
attenuation  rate  is  best  described  by  planar  divergence.  The  two-dimensional 
calculations  were  performed  using  the  finite  element  code  SAMSON  (Refer¬ 
ence  8).  SAMSON  is  based  on  a  small-strain  formulation  and  was  developed 
primarily  for  soil  structure  interaction  analysis.  The  finite  element  mesh 
for  this  problem  was  axisynmetric  about  the  centerline  of  the  cylindrical 
explosive.  Triangular  and  quadrilateral  elements  were  used  in  the  mesh.  The 
boundaries  of  the  mesh  were  located  at  a  distance  sufficiently  removed  from 
the  target  points  to  avoid  reflections  returning  to  the  region  of  interest. 

Due  to  the  high  pressures  involved  in  the  explosive  elements,  severe  mesh 
distortions  occurred.  In  an  effort  to  reduce  these  distortions  the  elements 
immediate  to  the  explosive  elements  were  treated  as  linearly  elastic  with  a 
somewhat  arbitrary  increased  stiffness.  This  adjustment  may  be  rationalized 
on  the  basis  that  the  mathematical  models  for  the  sand  and  concrete  do  not 
correctly  represent  the  behavior  at  high  strain  rates.  The  explosive  was 
modeled  using  the  Jones-Wi lkins-Lee  Equation  of  State  (Reference  9);  the  AFWL 
engineering  model  was  used  for  the  sand. 

The  SAMSON  code  results  are  shown  along  with  the  experimental  results  in 
Figures  25c  and  25d.  Figure  25c  reveals  that  the  blast  pressure  attenuation 
can  be  reasonably  computed  by  means  of  the  code.  This  figure  further  demon¬ 
strates  that  the  predicted  values  of  normal  stress,  the  experimental  values, 
and  those  computed  via  SAMSON  are  in  reasonably  good  agreement.  The  calcula¬ 
tions  indicated  in  Figure  25d  support  the  data  and  reinforce  the  conclusion-- 
fchat  of  near  planar  divergence  for  the  tests  with  a  burster  slab. 

As  mentioned  in  Section  IV,  the  data  in  the  a,X-plane  do  not  contain  any 
gravity  effects.  That  is,  the  Wi/3  term  in  the  denominator  of  X  is  actually 
(V^  •  G3)1/3  and  the  numerator  i s  •  G ;  thus  the  gravity  terms  theoreti¬ 
cally  cancel.  A  numerical  analysis  of  the  o  versus  X  data  could  not  detect 
any  statistic  revealing  any  gravitational  effects.  However,  if  normal  stress 


67 


is  plotted  versus  range  (dimensional),  then  the  influence  of  gravity,  as  it 
affects  explosive  yield,  is  clearly  shown.  Figure  26  is  such  a  plot  for  the 
burster-slab  data.  This  figure  was  constructed  from  the  least-squares  equa¬ 
tion  for  the  burster-slab  data  and  it  is  based  on  the  assumption  that  the 
phenomenon  under  study  is  described  by  the  surface  shown  in  Figure  6,  Sec¬ 
tion  III.  As  noted  in  Section  IV,  this  assumption  is  valid,  except  at  x  =  0.5 
due  to  an  oversight  in  the  experimental  design.  The  normal  stress  in  Fig¬ 
ure  26  is  defined  by 

a  =  9326  (R)-1*325  •  [G(W)1/3 j1  -241  (6) 

where 

o  =  normal  stress,  lb/ft2 
R  »  range,  feet 

G(W) 1 /3  =  gravity-scaled  yield,  ft-s-2  •  lb1/3 

Figure  27  is  a  plot  of  scaled  duration  and  impulse  from  selected  stress 
histories.  Only  a  few  of  the  data  traces  were  suitable  for  determination  of 
the  characteristics  (duration,  tD,  and  impulse,  i)  of  the  first  positive 
cycle.  Many  of  the  records  contained  reflected  impulses  and  most  exhibited  a 
zero  offset. 

The  slope  of  the  scaled  duration  relationship  is  in  essential  agreement 
with  near-field  data  presented  by  Higgins  et  al.  (Reference  3)  for  a  planar 
explosive  array.  The  increase  in  tp  in  the  near-field  (x  <  6)  is  apparently 
due  to  the  rapid  attenuation  of  the  high-frequency  components  of  the  stress 
wave.  In  the  far  field  (X  >  10)  the  stress  wave  should  propagate  at  a  near 
constant  frequency;  thus,  for  a  planar  configuration,  the  functional  relation¬ 
ship  between  scaled  time  (tQ  •  W“1/3)  and  scaled  distance  (x)  would  most 
likely  be  a  growth  curve  having  as  an  upper  bound  a  value  of  tQ  •  W-1/3  on  the 
order  of  0.1. 

The  following  equations  were  derived  from  Figure  27: 

t0  •  W1/3  =  0.0075  X°  -81+5  (7) 


T  .  U-l/3 


11  7  \  -0  -332 


where 


tg  •  W-i/3  =  scaled  duration  of  the  first  positive  cycle,  s-lb_1/3 

I  •  W-1/3  =  scaled  impulse  of  the  first  positive  cycle 
lb/in2 -s  •  lb-i/3 

x  =  scaled  distance,  ft-lb~1/3 

CONCLUSIONS 

The  purpose  of  this  study  was  to  evaluate  the  suitability  of  a  centrifuge 
as  a  device  on  which  blast  parameters  could  be  evaluated  on  a  small  scale. 

The  particular  problem  studied  was  concerned  with  the  simulation  of  blast 
parameters  in  a  sand  continuum  beneath  a  reinforced  concrete  burster  slab. 

Due  to  the  state  of  the  art  in  miniature  instrumentation  development,  the  only 
blast  parameter  that  could  be  evaluated  was  normal  stress.  Normal  stress  was 
measured  using  a  carbon-foil  pressure  sensor  mounted  on  stainless  steel  shim 
stock.  In  reality  the  development  of  a  normal  stress  gage  was  part  of  this 
study.  The  particular  gage  configuration  selected  was  based  primarily  upon 
size  considerations  to  be  consistent  with  the  concept  of  a  point  measurement 
in  a  high  G-level  acceleration  field.  Following  a  trial  development  period, 
during  which  the  final  gage  configuration  evolved,  the  gage  performed 
successfully. 

The  convincing  proof  of  the  suitability  of  the  centrifuge  as  a  blast- 
parameter  simulator  is  revealed  by  the  results  presented  in  Figure  25a.  Here 
the  data  for  the  tests  without  a  burster  slab  present  are  in  essential  agree¬ 
ment  with  the  preshot  predictions.  The  predictions  were  derived  from  an 
extensive  data  base  collected  over  several  decades  of  ground  motion  research. 

The  results  presented  for  the  tests,  which  included  a  burster  slab, 
reveal  a  signficantly  different  attenuation  rate  than  that  determined  for  the 
tests  without  a  burster  slab.  The  replication  in  the  data  does  not  support 
the  conclusion  that  this  difference  is  due  to  experimental  error  or  gravita¬ 
tional  effects.  Rather,  the  differences  are  due  to  the  burster  slab,  which 
alters  the  shape  of  the  wavelength  from  spherical  to  planar  divergence.  The 
phenomenological  differences  observed  are  further  supported  by  the  calcula- 
tional  effort,  which  indicated  significantly  different  attenuation  rates  for 
the  two  test  configurations. 


The  use  of  a  centrifuge  as  a  simulator  for  the  determination1  of  blast 
parameters  appears  to  be  a  workable  concept.  The  simulation  of  large  high- 
explosive  effects  through  gravity  scaling  permits  the  use  of  small  charges  in 
an  experimetal  device  that  is  easily  refurbished  after  each  test.  A  larger 
centrifuge  would  permit  the  use  of  larger  test  objects  and  reduce  the  effects 
of  reflected  waves  on  the  pressure-time  histories.  Also,  a  newer  centrifuge 
could  be  properly  instrumented  to  handle  low-voltage  data  channels  and  to 
isolate  them  from  deleterious  high-frequency  effects.  These  problems  notwith¬ 
standing,  this  study  does  demonstrate  that  a  centrifuge  can  be  used  to  evalu¬ 
ate  free-field  blast  parameters  (normal  stress). 
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Notify  ordnance  and  centrifuge  personnel  of  two 
impending  shots  on  centrifuge  in  Building  1001  on 
Kirtland  West,  tests  to  be  conducted  on  following 
day. 

Notify  explosive  safety  office  of  shots  in 
Building  1001  on  following  day. 

Take  documentary  photos  of  test  preparation. 

Transport  test  fixtures  from  CERF  main  facility  to 
Building  1001. 

Install  the  first  of  two  (2)  fixtures  on  centrifuge. 
Secure  all  instrumentation  cables  on  centrifuge  arms 
with  nylon  filament  tape  in  such  a  manner  that 
cables  will  not  stretch  and  break  under  strain  of 
high  (100  g)  forces. 

Hook  up  and  check  out  all  instrumentation  channels. 

Dry  run  system— electrical  and  explosive.  Post 
warning  signs— clear  area. 

Install  firing  module  on  centrifuge  and  place 
detonator  in  correct  test  position. 

Take  preshot  photos  of  completed,  ready- for-test 
configuration. 

Clear  out  centrifuge  area  of  personnel,  actuate 
safety  switch  on  centrifuge,  and  close  centrifuge 
pit.  Ensure  that  door  safety  switches  are 
actuated. 

Verify  that  all  systems  are  ready  for  test  sequence. 
Clear  area  topside  of  all  personnel  not  involved  in 
test. 

Energize  pumps  and  start  centrifuge.  When 
centrifuge  has  reached  correct  force,  fire 
detonator. 

Shut  down  centrifuge  and  turn  off  all  pumps.  Verify 
centrifuge  has  come  to  a  complete  stop. 

Verify  centrifuge  pit  area  is  free  from  explosive 
hazards  and  safe  for  personnel  to  enter.  Turn  off 
safety  switch  on  centrifuge. 
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T  + 

7  min 

Photo, 
EMCS,  TC 

Take  postshot  photographs,  inspect  test  fixture 
prior  to  disassembly  and  reassembly  for  second  test 
of  sequence. 

T  + 

30  min 

EMCS 

Begin  sequence  from  T  -  1  1/2-hour  step. 

T  +  30  min 

EMCS 

After  final  test,  notify  outside  personnel  testing 
is  completed  for  this  day. 

The  word  HOLD  from  any  test  personnel  will  cause  the  countdown  to  be 
stopped.  Test  will  resume  from  T  -  5  min  countdown. 

Hold  conditions:  Any  malfunction  of  either  the  centrifuge, 
electronic,  or  explosive  systems  which 
would  render  the  test  procedure  unsafe,  or 
result  in  unacceptable  test  conditions 

Responsible  Officers:  TC  -  J.  Nielsen 

EMCS  -  S.  Scales 
CENT  -  P.  Adams 
ELEC  -  G.  Stewart 
PHOTO  -  M.  Davila 
P.I.  -  J.  Nielsen 
ORD  -  R.  Venegas 

Prepared  By:  _  Approved  By:  _ 

Date:  Date: 
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